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Abstract 



We discuss the nonholonomic Chaplygin and the Borisov-Mamaev-Fedorov systems, for 
which symplectic forms are different deformations of the square root from the corresponding 
invariant volume form. In both cases second Poisson bivectors are determined by L- 
tensors with non-zero torsion on the configurational space, in contrast with the well known 
Eisenhart-Benenti and Turiel constructions. 

1 Introduction 

At the beginning of the 20th century SA. Chaplygin showed that a two degree of freedom non- 
holonomic system possessing invariant measure can be reduced to Hamiltonian after a suitable 
change of time [TU1 [H] - Using this process of the Chaplygin hamiltonization we can get the usual 
Hamilton- Jacobi equation, variables of separation, the Abel-Jacobi equations, quadratures, etc 



In bi-hamiltonian geometry separability is invariant geometric property of the distribution 
defined by mutually commuting independent functions Hi, ... , H n [201 H21 [23] . In fact there is 
neither Hamilton- Jacobi equation, no time which describes only some partial parametrization 
of geometric objects. So, in this paper we want to show how these standard bi-Hamiltonian 
geometric methods may be directly applied to nonholonomic systems without any change of 
time. The second aim is to discuss a deformation of the Turiel construction [55], which ap- 
pears only in the nonholonomic case and gives rise to some interesting modifications of the 
standard conformal Killing tensors that lie at the heart of classical Eisenhart-Benenti theory of 
separability 0[3J[TI]. 

As an example, we will consider a rolling of dynamically asymmetric and balanced ball 
over an absolutely rough fixed sphere with radius a. At a — > oo one gets a Chaplygin problem 
on a non- homogeneous sphere rolling over a horizontal plane without slipping [10J. Thus, we 
are able to compare bi-Hamiltonian methods with the the Chaplygin hamiltonization. 

Let u — (uii,u)2,u>3) be an angular velocity vector of the rolling ball. Its mass, inertia 
tensor and radius will be denoted by to, I = diag(7i, I2, h) and b, respectively. According to 
[5J, the angular momentum M — (Mi, M2, M3) of the ball with respect to the contact point 
with the sphere is equal to 



Here 7 = (71, 72, 73) is the unit normal vector to the fixed sphere at the contact point, E is the 
unit matrix and (., .) means the standard scalar product in R 3 . All these vectors are expressed 
in the so-called body frame, which is firmly attached to the ball, and its axes coincide with the 
principal inertia axes of the ball. 

For further use we rewrite the relation (jl.ip in the following equivalent form 



m HI E MM HZI- 



M = (I + <2E) w 



d(7,w)7, 



d = mb 2 . 



(1.1) 




(1.2) 



and 



where 

Oi 
A= ( a 2 ) =(I + dE)- 1 , 
a 3 

According to [TUJ [5] , there is no slip nonholonomic constraint associated with the zero velocity 
in the point of contact. It allows us to reduce equation of motion to the following form 

M = M x uj , 7 = K7 x lo . (1.4) 

where x means vector product in R 3 and k — a /(a + b). 
At any k there are three integrals of motion 

H x = {M,u), H 2 = (M,M), Ci = (7,7), (1-5) 

and invariant measure 

^=^jdjdM. (1.6) 
If k = ±1 there is one more integral of motion 



bi 






















b 3 



C 2 = ( 7 ,BM), B= 6 2 ^trA^ + C/c-lJA- 1 . (1.7) 

V b 3 J 

At d = we have 3(7) = and a; = AM. In this case \x is a standard volume form with 
constant density and, therefore, equations (11.41) describe some hamiltonian flow. Namely, at 
k = 1 equations (jl.4l) can be identified with the Euler-Poisson equations describing the rotation 
of a rigid body around a fixed point, whereas at k = — 1 equations (|1.4[) describe the so-called 
Contensou model of Fleuriais gyroscope [J2] . 

At d equations (|1.4p describe the rolling of a dynamically nonsymmetric sphere over a 
horizontal plane or a fixed sphere without slipping. The case k = 1 is the so-called Chaplygin 
system [TU] and case n = —1 will be referred as the Borisov-Mamaev-Fedorov system. A detailed 
description of these nonholonomic systems may be found in [5J [6j [7l [HJ [9l [13] . 

In pj |H1 [H] authors change time variable in the equations of motion (|1.4[) first and only 
then study the Poisson structure of the resulting equations. In the Hamiltonian mechanics 
transformation of time can drastically change almost all the invariant geometric properties of 
the initial system, such as the Lagrangian foliation, compatible Poisson structures, Lax and 
r-matrices, bi-Hamiltonian construction of the variables of separation, etc [TU [19]. For the 
nonholonomic system it can change even the initial Hamilton function [TJ 1131 117] . 

Our main aim is to get a family of Poisson brackets associated with the invariant measure 
(|1.6p in framework of the bi-Hamiltonian geometry, i.e. without any change of time. It means 
that we for a while forget about the equations of motion (|1.4I) and try to solve the following 
geometric equations 

Pdd, 2 = 0, (PdH u dH 2 ) = {H 1 ,H 2 } = 0, [P,P]=0, (1.8) 

where [., .] is the Schouten bracket, with respect to the Poisson bivector P. 

So, in our approach the Hamiltonization process is equivalent to a search of the Poisson 
structure satisfying to equations (|1.8[) . i.e. 

Hamiltonization ■<=>■ Poisson bracket 

because using this bracket we can always get new Hamiltonian system 

-^pz k = {Hi,z k } , 

with new time t' and the same integrals of motion in involution. 
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Remark 1 In our case at d = we have Hamiltonian system with canonical Poisson bracket. 
So, Hamiltonization is equivalent to existence of the proper deformations of this canonical 
Poisson bracket. Obstacles to such deformations are well known, see, for instance, the geometric 
quantization theory. 

The principal disadvantage is that equations (|1.8[) have infinitely many solutions [20, 21, 22 . 
So, in order to get any particular solution we have to set aside an invariance and to narrow the 
search space. 

1.1 Spherical coordinates 

In order to to narrow the search space we will use the notion of natural Poisson bivectors on 
the Riemannian manifolds [22 . In this case we have to reduce our initial phase space to the 
cotangent bundle of the unit two-dimensional Poisson sphere. 

Namely, we can avoid the solution of the first equations PdC\ t 2 = in (|1.8[) using the 
slightly modified Euler variables 




such as 

C 1 = ( 7 , 7 ) = l, C a = ( 7> BM)=p^ 

In the generic case the remaining equations in (|1.8[) have infinitely many solutions [20 . In order 
to find at list one particular solution we impose the following restriction 

C 2 = ( 7 ,BM) = Pi> -0. (1.10) 

In this case we have natural integrals of motion (|1.5j) and, therefore, we can solve our 

equations using the notion of natural Poisson bivectors [22] . 

At K = 1 coordinates ((f), 9) in (|1.9[) are usual spherical coordinates on the unit sphere S 2 

at 

K = 1 , b\ = b 2 = b 3 = tr A -1 = 1 . 

At k = — 1 we replace parameters ai and J; = a^ 1 on 6, 

2 2 2 

K = ~ 1 ' °i — 7 — rTT' ° 2 = 1 — TIT' «3 = 7 — —r-, 

02 + o 3 bi + b 3 oi + o 2 

in order to get more foreseeable formulas. 

After the change of variables (|1.9I) . at n = ±1 we obtain two different dynamical systems 
on the common phase space M. which is topologically equivalent to cotangent bundle T*S 2 
to the sphere. These systems have a common invariant volume form (|1.6p and two different 
Poisson structures, see next Section. 

In order to show the difference between the Chaplygin and Borisov-Mamaev-Fedorov sys- 
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terns we present one equation of motion 

• (m — 0,2) sin20 (cosO g(<i3 — a\ sin 2 — 02 cos 2 0) sin 20 



2 V sin ^ 
g(ai — ei2) 2 sin 2 20 sin 2 6* 



ai cos 



02 sm 



k = -1 



1 



63(61 — 62) sin 20 cos # 



sm f 



M>2(&l + &3)(&2 + & 3 ) 

g(6i - 6 2 )(6 3 (6 2 cos 2 + 6| sin 2 0) - 6 2 6?.) sin 20 sin 26*' 



63(61 +h)(bi + 6 3 )(6 2 + 63) 
2(63(61 sin 2 + 6 2 cos 2 0) + 6i6 2 ) - 
and one integral of motion 



g6 3 (6i - 6 2 ) 2 



(61 + 63) (6 2 + 6 3 ) sin 2 20 sin 2 



K = 1 , 



H2 = -rjTEPl+Pei 



sm 



Ho = 



(bj cos 2 (f> + bl sin 2 0) cos 2 6> 1 



6?6t sin 2 ( 



+ 6?W + 



(6 2 — 63) sin 20 cos ( 
b\b\ sin# 



Po 



■P<pPe 



6 2 sin 2 + 6| cos 2 2 
+ 6f6l 

Of course, any calculations for the Borisov-Mamaev-Fedorov systems require more efforts and 
large-scale resources in comparison to the same calculations for the Chaplygin system. 



2 Invariant measure and Poisson brackets 



Let M. be a smooth symplectic manifold endowed with a symplectic form fl which in the 
Darboux coordinates 

z = (q,p) = (qi,...,q n ,pi,...,p n ) 

reads as 

fi = dpi A dqi + . . . dp n A dq n . (2-1) 

The volume form Q 2 on Ai is invariant under all hamiltonian diffeomorphisms by the Liouville 
theorem. 

If we have another invariant volume form \i on the same manifold M. , we can get another 
symplectic form taking a formal square root on /x, because 



A* = 



(2.2) 



However, in our case invariant volume form [i = ^/g A (|1.6p is invariant with respect to the non- 
hamiltonian flow (|1.4|) and, therefore, we have to deform its formal square root (|2.2|) . We will 
describe these deformations using Poisson bivectors, instead of the corresponding symplectic 
forms. 

We rewrite the Poisson bivector P associated with the canonical symplectic form SI (|2.1I) 
in the following tensor form 







P = 



E 



dL ki dL 



\ 



dqj dqi 



Pk 




(2.3) 
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where L is an identity (1,1) tensor field on a configurational space. 

We use such unusual notation because any torsionless (1,1) tensor field L'(qi, . . . , q n ) on a 
configurational space Q with coordinates qx,...,q n determines another Poisson bivector 



P' 







\ 



fdLL dL> 



\ 



\ dqj dqt 



Pk 



fe=i 

on M. , according to [35] . The corresponding Poisson brackets read as 
= 0, {QuPjY = L'ij > 



(2.4) 



ki 



dL 



kj 



dqi 



Pk ■ 



The vanishing of L' torsion entails that P' (12 .4[) is a Poisson bivector compatible with P, i.e. 

[P,P] = [P',P] = [P',P'} = 0. 
The torsion of the (1,1) tensor field A equals to zero, if for any vector fields X, Y 



T A {X, Y) = C AX AY - a(l ax Y + C A yX- A£ x Y 
Here Cx means the Lie derivative along X. 



0, 



Remark 2 In the framework of the Eisenhart-Benenti theory tensor field V yields special 
conformal Killing tensor called the Benenti L-tensor, Killing-Stackel space, Stackel web etc 
[3J|3]. Of course, we can try to transfer the corresponding geometric machinery to nonholonomic 
theory. 

One possible generalization of the Turiel construction (|2.4[) has been proposed in [22 . Here 
we consider some other generalizations related with the nonholonomic integrable systems. 

In our case n = 2 and the Darboux coordinates on A4 = T*S 2 are standard spherical 
coordinates (|1.9[) , thus, 



Ql = 4> , Q2 = , pi = p,p , 



P2 



(2.5) 



At d = and k = ±1 we have the constant invariant measure and, therefore, integrals of motion 
-Hi, 2 (|l-5p are in involution 

{H u H 2 } = 0, d = 0, k = ±1, 

with respect to the Poisson brackets associated with the canonical bivector P 

2.1 Case k = 1 

At d and k = 1 substituting another torsionless tensor field 



L a = 



1 



1 



1 
1 



(2.6) 



into the definitions (|2.4I) , one gets the desired solution of the equations ([1.8 



P 9 = 



( 1 

* 

* * 



1 /31ng 

2 { 86 1 



ding 



-po 



(2.7) 
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Remark 3 This deformation has a similar form with the well-known relation between the 
modular vector fields 

Xgfi =Xfj,- X\ ug 

associated with the volume forms /i and v = g/j, |15[ 126] . 

Remark 4 At d = function g (|1 .3[) equals to unit and, therefore, at this limit one gets 
standard canonical Poisson bivector P — lim P a . 

In terms of initial variables (7, M) this Poisson bivector P g (|2.7p has been obtained in [5]. 

Proposition 1 JfJj] Integrals of motion if 1.2 HI .5]) are in involution with respect to the Poisson 
bracket associated with the Poisson bivector P g {2.7}) 

{H u H 2 } g = 0, d>0, k = 1. 

The corresponding volume form 

v = Pg 2 = -2gdgdp 

is invariant with respect to a Hamiltonian flow associated with new time t g defined by 

— z k = {Hi, z k } g , fc = l,...,4. 
atg 

We can easily relate new and old time variables 

At g ~ Vgdi (2.8) 
because at k = 1 initial equations of motion are equal to 

± Zk = ^{ Hl ,z k } g . (2.9) 

Transformation of time (|2.8[) has been proposed by Chaplygin in jlOJ . Namely this process is 
to be referred to as the Chaplygin Hamiltonization, see IH [9j [13l [17] . 

Remark 5 One of the global invariants in Poisson geometry is a modular class. It is an 
obstruction to the existence of a measure in Ai which is invariant under all hamiltonian flows 
[13 [IHl ■ So, in fact it is a geometric obstruction to the Hamiltonization process. 

For the manifold M. endowed with a Poisson bivector P, its modular class is an element of 
the first Poisson cohomology group. In Section 3 we discuss some elements of the second Poisson 
cohomology group and the corresponding Poisson bivectors P' compatible with P, which allows 
us to get variables of separation without Hamiltonization. 

2.2 Case k = -1 

It is easy to see, that at k = — 1 the integrals of motion H1.2 (|1.5p do not commute with respect 
to the Poisson brackets associated with bivector P g ()2.7[) 

{H 1 ,H 2 } g ^0, d>0, k = -1. 

So, we have to propose another deformation of the canonical Poisson structure applicable to 
the Borisov-Mamaev-Fedorov system. 

Let us try to solve our geometric equations 

(PdH 1 ,dH 2 ) = {H 1 ,H 2 } = 0, [P,P}=0, (2.10) 

by "brute force" method, using similar to ()2.7[) anzats 

/ f(cf>,d) \ 

p= * h(cf>,6) 

* * u((f>,9)p$ + v((j>,9)po 

\ * * * J 

As a result we have the following 
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Proposition 2 At k = — 1 the integrals of motion Hx,2 U.5\) are in involution with respect to 
the Poisson brackets associated with the Poisson bivector 







L, 



\ 



fc=i 



dL 



nkj 



(1 + ??) % 



Pk 



(2.11) 



where 



L, 




and r\ 



2sm 2 e(bj - (&! + 62)63 + hb 2 
b 2 3 (d~ 1 {bi + b 2 )-2) 



(2.12) 



Here L v is the (1,1) tensor field with non-zero torsion, in contrast with the tensor field from 
the Turiel construction {2.$. 



The proof is straightforward. 

Tensor field L v may be considered as an additional deformation L g (|2 .6[) by function 77 
depending only on variable 6, parameters d and bk, such as 



lim n = 



lim P n = P. 



Moreover, -q = for the axially symmetric ball at 63 = b\ or 63 = 62- 

Remark 6 At present we don't know any physical meaning of the function r)(6) and the geo- 
metric explanation of the deformation (|2.11[) . It will be interesting to understand the relations 
between L v and the theory of Killing tensors with non-zero torsion. 

The Poisson bivector P n (|2.11l) may be rewritten as follows 



/ 1 

* 



VI 



* * 



V 



(1 + v) 





(1 + v) 

ding d In g 



* * * 



The corresponding volume form 



P' 



2g 



■ dqdp 



is a more complicated deformation of the invariant volume form (jl.6[) introduced in |27j 

/i = Vg dgdp . 

This new volume form is invariant with respect to a Hamiltonian evolution associated with new 
time trj defined by 

z k = {Hi, z k } v , k = l, ...,4. 

Relation between the initial and new time variables is also more complicated then in 
because at k = — 1 initial equations of motion (|1 -4[) read as 



At Zk = Vv 1 + &2 + 63 + Wl J^ jffl ' Zfe ^ _ Vg[l + w 2 ){H 2 ,z k } n . (2.13) 
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Here 

rj (h + b 2 ) {bs{bi cos 2 4> + b 2 sin 2 4>) - &i& 2 ) 



i/'i 



U?2 



l + V {bi - h)(b 2 - 63) 



(2.14) 



7/ ^3(61 cos 2 + 62 sin 2 <fi) — b\b 



'2 



l + V {bi-b 3 ){b 2 



So, without the intermediate time transformation we did not get conformally Hamiltonian 
system from 0|5], because we consider, in fact, two different systems with geometric point of 
view, see [4] 1131 117] . The modern theory of conformally Hamiltonian systems may be found in 



3 Second Poisson brackets 

In this Section we want to get another solution P' of the equations (|1.8|2.10p . which is com- 
patible with the first solution P obtained earlier, i.e. 

[P,P']=0. 

Compatible bivectors P' are the 2-cocycles in the Poisson cohomology defined by P on the 
Poisson manifold M. , whereas the Lie derivatives of P along vector field X 

P' = C X P 

are2-coboundaries. So, in order to get the desired solution of (|1.8I2.10[) we will use the Lie 
derivatives along the vector fields X with linear in momenta entries. 

In bi-Hamiltonian geometry equations of motion usually have the following form 

— z k = si {H u z k }' + s 2 {H 2 ,z k }' , (3.1) 
at 

where {., .}' is the second Poisson bracket associated with P' and S1.2 are some functions on 
dynamical variables [35] . 

If s± — and s 2 = const we have a bi-Hamiltonian dynamical system. If si = and s 2 is 
arbitrary, one gets the so-called quasi bi-Hamiltonian system. At s\. 2 7^ we have bi-integrable 
dynamical system (20] [22]. So, the equations of motion (|2.13|) for the Borisov-Mamaev-Fedorov 
system have the standard bi-Hamiltonian form. 



3.1 Case d = and k 



At d — 1 we have the hamiltonian flow (|1.4p associated with the canonical Poisson bivector P 
(|2.3p . It is easy to prove that the integrals of motion H\ j2 (|l-5p are in bi-involution 



{H 1 ,H 2 } = {H 1 ,H 2 }' = 0. 



0, 



with respect to canonical Poisson brackets associated with bivectors P (|2.3p and P' (|2.4p de- 
termined by the following (1,1) torsionless tensor field 22 : 



11 = 



( 



\ 



a\ cos (f> + a 2 sm 
(cii — a 2 ) sin 2(j) 



(ai — a 2 ) sin 20 cos( 



s,m6 



cos 6 sin 9 sin 9 + (ai sin <j) + a 2 cos <j)) cos 9 



(3.2) 



The Turiel bivector P' (|2.4p may be rewritten as the Lie derivative P' — Cy P of the canonical 
bivector P along the vector field Y = Y^dj with the following entries 

Y*)-- L [po 



Y 1 ' 2 =0, 



(3.3) 
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Here V stands for the transpose of the matrix V . 

Remark 7 The corresponding volume form A' = P' 2 is invariant with respect to the new 
time defined by 

-^ i Zk = {Hi,z k }' , fc = l,...,4. 

It is neither bi-Hamiltonian nor quasi bi-Hamiltonian system |22] and, therefore, this repara- 
metrization of time looks like the Hamiltonization for the Borisov-Fcdorov system. 

The eigenvalues u,v of the recursion operator TV = P'P^ 1 are the roots of the following 
polynomial 

B(X) = (A - u)(X -v) = X 2 -tr (L'L- 1 ) A + ^j-^ = . (3.4) 

Of course, in this case coordinates u, v are the standard elliptic coordinates on the sphere 
defined by 



(X-u)(X-v) 7? 7| t| 



(A - ai)(A - a 2 )(A - a 3 ) A - ai A - a 2 A - a 3 



(3.5) 



3.2 Case d = and k = — 1 

At g? = and k = — 1 the integrals of motion H± t 2 (|1.5p are in bi-involution with respect to 
canonical Poisson bracket and the second bracket associated with bivector P' (|2.4[) defined by 
the following (1,1) tensor field 



L' = 



V 



ci cos <p + C2 sm 
(ci — C2) sin 20 



(ci — C2) sin 20 cos 9 



cos sin 8 C3 sin + (ci sin 



sin 



C2 cos 2 (/>) cos 2 1 



/ 



(3.6) 



where Cj = di/bi. The eigenvalues of the recursion operator iV = P'P 1 coincide with standard 
elliptic coordinates on the sphere 



7t 



(U - Cj)(v - Ci) 
(Cj - Ci)(c k - Ci) 



(3.7) 



The conjugated momenta p u ,Pv are defined by standard relations 

1 2e i jkJj'Yk(cj - c k ) 



Mi 



■Ufii - u)p u - {ci ~ v)pSj 



(3.8) 



where Sij m is a completely antisymmetric tensor. 

In terms of these Darboux-Nijenhuis variables u,v and p ui p v , our Poisson bivectors look 

like 



P = 



/ 1 \ 
1 
-10 

\ -100/ 



P' = 



( u Nj 

v 

-u 

\ -v / 



whereas in terms of initial physical variables, first bivector P reads as 



/ 

* 



P = 



s(S * * 
* * * 





"&2&373 



* * b 3 b 2 -f2 

* * * 



M373 -^1^272 \ 

626171 
-636171 



6§M 3 




-6 2 M 2 
b\M x 



(3.9) 
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It is evident that any functions fi(u) and f 2 (v) are variables of separation as well. So, the 
following trivial point transformation 



(3.10) 



preserves the separability property of distribution defined by the functions if 1,2- Namely, 
according to [20] , these integrals are in involution with respect to the Poisson brackets associated 
with the following bivectors 



P' = 



V 









h{u) 














f 2 (v 
















-h(v) 









/ 



(3.11) 



Of course, in terms of initial physical variables these bivectors have more complicated form. 
For instance, tensor field 

' P 
p 



L" 



(3.12) 



where 



c 



cos 2 9 



63(62 + hi) cos 4 9 



(6i6 2 + 63(61 sin 2 + 6 2 cos 2 </>)) 
a „„ 2 a _l 



cos 2 6 sin 2 1 



sin 2 6 cos 2 1 



6162 + 63(61 sin 2 c\> + b 2 cos 2 0) 61(62 + 62) 62(61+63) 

yields bivector (|3.11[) associated with new variables of separation (|3.10[) defined by 

2( u 6x(6 2 + 63) - 2) (t*6a(6i + 63) - 2) 



A(«) = 

f 2 (v) = 



w6i6 2 (6i + 6 3 )(6 2 + 6 3 )(u6 3 (6i + 6 2 ) - 2) 

2(^1(62 + 63) -2)^62(61+63) -2) 
vbib 2 {bi + 6 3 )(6 2 + 6 3 )(w6 3 (6i +6 2 ) - 2) 



It is natural that the canonical transformations (|3.10p preserve the first bivector and change 
the second bivector P' simultaneously with coefficients S\ 2 in the equations of motion (|3.1[) . 
Of course, some geometric properties of these equations are invariant with respect to such 
transformations. 

Proposition 3 At d = and k = —1 there does not exist nontrivial linear in momenta Poisson 
bivector P" , which is compatible with the canonical ones, such that s± = in \3. 1}) . 

So, at k = —1 the dynamical system \l-4\ ) is only bi-integrable, whereas at n = 1 it is 
bi-Hamiltonian. 

By adding equations (|3 . 1 1) with s± = and compatibility condition [P, P"] — to the initial 
equations (|1.8|2.10|) one gets an overdetermined system of algebro-differential equations. If the 
entries of P" are linear nonhomogeneous polynomials in momenta, then the system has only 
trivial solution P" = 0. 

Remark 8 According to [5], at d = dynamical systems with k = ±1 are related to each other 
by the Poisson map M —> BM and the trivial change of time 

£-> -t . 

In Proposition 3 we proved that even such seemingly harmless transformation leads to a loss of 
very important geometric property. Namely, after this change of time the new system becomes 
non bi-Hamiltonian with respect to initial integrals of motion. 
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3.3 Chaplygin system, k = 1 

According to [23], let us introduce the vector field X = ^X^dj with the following entries 



X 1 = 0, 



x l+3 = 



7 x A g ( 7 x M) 



i = 1,2,3. 



(3.13) 



where A„ is the following 3x3 matrix 



A g = A + dg(j) A (7 (g) 7) A . 



entering into the angular velocity definition (ll.2[) . 



Proposition 4 f^f TTie Lie derivative of P g \2. 7| ) along the vector field X i3.13\) is the desired 
second solution of the equations U.8l2.l0\) compatible with the first solution 

Pg =£xP g , (3.14) 

so that the integrals of motion Hip \l-5\) are in bi-involution 

{Hi, H 2 } g = {Hi, HzYg = , (3.15) 

with respect to a pair of the corresponding compatible Poisson brackets. 

In spherical coordinates this bivector looks like a deformation of the Turiel construction (|2.4[) 



/ 



P' 

g 







L' 



\ 



\ 



-L' 



3ij 



dL 9k 3 \ 

Vkj Q Pk 



(3.16) 



fc=i 



Similar to the tensor field L v (|2. 12[) , this (1,1) tensor field 
d y/g sin 2 



L' 



fgL' 



1 — daa 



0102-03(^1 cos 0+O2 sin </>) 

— 03+03(01 sin 2 (/>+02 cos 2 



has a non-zero torsion too. 

Functions Xki , Ukj depending only on the coordinates <f>, 9 can be easily restored from the 
relation (|3.14l) . which in spherical coordinates reads as 



P'=C z P g , Z^ 2 = 0, 



Z 3 
Z 4 



p,p 
Po 



(3.17) 



Here L' g stands for the transpose of the matrix L' g . 

Using this second Poisson structure for the Chaplygin system we can rewrite the equations 
of motion (|1.4|) in the following form 



where 



df 



si = -1 - da 3 + 



d ( 

Zk = - 



2 — — {#l,Zfc} a + {.ff2,*fe} fl J , 



as sin 



(03(01 + 02) — 0102) cos 2 6 — d(da\a2 — di — 02)03 



dai0203 — 0102 cos 2 — 03(01 cos 2 <fi + 02 sin 0) sin 9 



The eigenvalues it s , u g of the recursion operator N g = P' g P g 1 are defined by the relation 



(X-Ug)(X-Vg) 

(A - oi)(A - a 2 )(A - a 3 ) 



5(7) 



7 2 (l-doi) 7 2 2 (l-^2) , 7 3 2 (l-^3) 



A — ai 



A - a 2 



A - o 3 



(3.18) 
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which bears a resemblance to the usual definition p.5[) of the elliptic coordinates on the sphere. 
Namely these variables have been obtained by Chaplygin who used hamiltonization [TU] . 

As above, using the trivial point transformations (|3.10p we can change this second Poisson 
bracket and the equations of motion (|3.1| associated with this bracket. For instance, let us 
consider another (1,1) tensor field 



jj> — _J_ 



L' 



1 



1 — da 3 



where 



da\a 2 a 3 + a\(i2 cos 2 8 + 0,3(0,1 cos 2 





Ml 




Pi ) 


a 2 sin 2 


<j)) sin 2 



(3.19) 



Pi 



= I da* 



sin 2 <h cos 2 ( 



cos' 



+ ( 



da? 



cos sm 



1 ) a 2 — a 3 sin 



Pi 



pi + dsin 2 6(a\ — a 3 )(a 2 — a 3 ) . 



Substituting this tensor field into the Lie derivative p. 171) , one gets a new Poisson bivector P' g ' 
with the following properties. 

Proposition 5 At k = 1 i/ie initial equations of motion \1.J$ have the following ft 



orm 



d 7 _ Vlrrr rl _ (1 ~ da 3 )g r „ y, 



(3.20) 



ociated with bivector P' 1 . 



dt " 2 

where {•> -}g *s Poisson bracket a.ss 

So, equations of motion for the nonholonomic Chaplygin system are conformally Hamiltonian 
equations with respect to both the first bracket {., .} g with first integral of motion Hi (|2.9| and 
the second bracket {., .}' g ' with second integral of motion H 2 (|3.20p . 



3.4 Borisov-Mamaev-Fedorov system, « = — 1 

As usual, at k = — 1 there exist many linear in momenta solutions of the equations (|2.10l) , which 
are related to each other by point canonical transformations — > /i(A a ) p.lOp . where Aj are 
the Darboux-Nijenhuis coordinates, i.e. the eigenvalues of the recursion operator. 
All these solutions have the form 



P' 







V 



-LL 



E 



dL 1 dL'. 

Vki Vkj 



dqi 



Pk 



(3.21) 



Let us consider only one solution associated with a relatively simple tensor field 



LL 



1 

vi 



L' 



(l + r))a + 





Mbi + h + h) 



(3.22) 



{bi+b 2 )(b 1 +b 3 ){b 2 + b 3 ) 

depending on the arbitrary number a. As above, this tensor field has a non zero torsion at 
generic a. 

Functions Xki,y k j depending only on the coordinates <j), 9 can be easily restored from the 
other definition of the same bivector 

p: = c z p n , 

where the entries of the vector field Z — ^2 dj are equal to 



Z 1 ' 2 = o. 



z 4 



1 

(l + n)- 1 



Pel, 

po 



(3.23) 
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Remark 9 It is easy to see, that at k = ±1 entries of the Liouville vector field Z can be 
rewritten in the common form 





where Z = (Z 1 , Z 2 , Z 3 , Z 4 ) is the vector of the entries vector field, whereas z — (qi, q 2 ,pi,p 2 ) 
is the vector of Darboux coordinates . Geometric origin of this new construction is unclear of 
yet. 

So, for the Borisov-Mamaev-Fedorov system we have a pair of compatible Poisson bivector 
on M and, therefore, this manifold is bi-Hamiltonian. 

Proposition 6 At n = —1 the integrals of motion H\ 2 hi. 5)) are in bi-involution 

{H 1 ,H 2 } v = {H 1 ,H 2 y ri = 0, (3.24) 
with respect to the Poisson brackets associated with bivectors P n and P' . 

Using the second Poisson brackets {., .}^, we can rewrite the initial equation of motion (jl.4[) in 
the standard form (I3.1[) with relatively big coefficients s% >2 . 

At a = in (|3.22[) . the eigenvalues and of the recursion operator N v 
defined by the relation 



(X - u v )(X - Vr,) 1 / 7l 2 7I /37| 



(A- ci)(A- c 2 )(A- c 3 - 5) C(t) \A-ci A-c 2 \-c 3 -6 
depending on two constants 

2d(h - 63) (62 - b 3 ) 
b 3 (b 3 -2d)(b 1 +b 2 ) + 2db 1 b 2 ' 

Ad hb 2 (h-h)(b 2 -b 3 ) 



(3.25) 



1 



(3.26) 



b 3 {b 3 - 2d)(b 1 + 6 2 ) + 2dbxb 2 b 3 (bx + b 2 )(h + b 3 ){b 2 + b 3 ) ' 

and one function on 9 = arccos73 

C(7 ) = 1 2d(6 1 -ft 3 )(6 a -ft 3 ) 2 

^" b 3 (b 3 - 2d)(h + b 2 ) + 2dhb 2 13 ' 

This relation is very close to (|3.5p and (j3. 18[) and, of course, at d = coincides with the 
definition (|3.7|) of the elliptic coordinates on the sphere § which are variables of separation for 
the corresponding Hamilton- Jacobi equation. 

Now we have to compare this Darboux- Nijenhuis coordinates with the variables of separa- 
tion obtained in [7J H] by hamiltonization process. 

Proposition 7 Variables of separation qi^ from are related with the eigenvalues u n and 

v n \3.25)) of the recursion operator H n by trivial point transformation similar to \3.10\) . 

Let us reproduce the definition of variables of separation qi t 2 from [S], see formulae (3.2): 



det I / (qi - a) (q 2 - c*) _ ^ 



y (J l - d)Jj J fc G(qi,q 2 ) y {cj - Ci){c k - Ci) 
where 

, . (bi + b 2 - 2d){h + 63 - 2d) (6 2 + b 3 - 2d) 

lqi ' q2j (bi+b 2 )(h+b 3 )(b 2 + b 3 ) 
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Substituting this definition into the equation (|3.25[) one gets the desired point transformation 



Ur) = F(q{) , v n = F(q 2 ) , 

where 

(b% + (bi + b 2 - 2d)b 2 3 + bib 2 b 3 + 2d6i6 2 )q - 4d 

= (b 1 +b 3 )(b 2 + b 3 )(d(b 1 b 2 q-2) + b 3 ) " 

Associated with the variables qi,2 tensor field L' in (|3.21[) is more complicated then tensor field 
(|3.22j) and, for brevity, we omit this expression. 



4 Separation of variables 



In geometry, instead of an additive separation of variables in the partial differential equation 
called the Hamilton- Jacobi equation, we have some invariant geometric property of the La- 
grangian distribution defined by n independent functions Hi , . . . , H n . 

Namely, an n-tuple Hi , . . . , H n of functionally independent functions defines a separable 
foliation on M, dinxM = n, if there are variables of separation (qi, . . . , q„, pi, . . . , p„) and n 
separated relations of the form 



<&i(qi, p,, #!,...,#„ 



, 



1, 



with det 



dHi 



(4.28) 



It simple means, the common level surfaces of Hi, ■ ■ ■ , H n form foliation and every leaf of this 
foliation may be represented as a direct product of one-dimensional geometric objects defined 
by separated relations (|4.28[) . Usually we have a direct product of n algebraic curves, because 
$i are polynomials in q^ and p,. 

It can be easily shown 20 , that condition (|4.28j) entails the involutivity of Hi with respect 
to the compatible Poisson brackets 



hi^j}f = {Pi> Pi}/ = °> {Pu q 3 }f = % /i(Pi. Qi) . 



(4.29) 



depending on the arbitrary functions fi, . ■ ■ , f n . In fact, this definition of separability implicitly 
appeared in the Lagrange proof of the Jacobi theorem, but both Lagrange and Jacobi used only 
canonical brackets 



{qt,qj} = {pi,pj} = o, {<u,Pj} = 5ij, (4.30) 

which belongs to the family (|4.29[) . 

In bi-Hamiltonian geometry eigenvalues qi , . . . , q n of the recursion operator are the desired 
coordinates of separation or the Darboux-Nijenhuis coordinates. It is a sequence of the fact 
that the distribution tangent to the foliation defined by Hi , . . . , H n is Lagrangian with respect 
to the symplectic form P _1 and invariant with respect to recursion operator N = P'P -1 . So, if 
we know these coordinates, then we have to explicitly find the conjugated momenta pi, . . . , p„ 
and the separated relations (|4.28[) . 

In the Chaplygin hamiltonization method momenta are defined by complete integrals S of 
the Hamilton- Jacobi equation after the corresponding change of time 

d 

Pj = Q^Sj(qj,ai,... f a n ). (4.31) 

On this step one usually gets momenta p^ which have more complicated Poisson brackets with 
coordinates q^ (|4.29|) instead of standard canonical brackets ([4.301) . 
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For instance, let us consider definition of the separated momenta from [9], see formulae 
(3.14): 



M = (J t - d) 2 y/lcj - qi)(cj - q£)y/(c k - qi)(c fc - q 2 ) 
Jfbi ' 2^G(qi,q 2 )( U - W ) 



P2 P2 



(qi - c,-)(qi - Cfc) (q 2 - Cj-)(q 2 - c k ) 

It is easy to observe, that at d = variables qi. 2 (13.27[) coincide with the standard elliptic 
coordinates u, v (|3.7p on the sphere S, but the corresponding momenta 

pi = 4>{u)p u , p 2 = 4>{v)p v (4.33) 

differ from the standard canonical variables p u and p v (13.81) on T*S. 

Moreover, after substituting 7$ p.27[) and Mj ()4.32[) into C 2 = b^Mi one gets C 2 7^ 
even at d = 0. So, we suppose that the definition of momenta in [S] contains some misprint 
and, therefore, we have to define these variables correctly. 

4.1 Chaplygin system, k = 1 

According to [23 [21] > we can use the following recurrence chain 

0i = {u, H k } g , (t>2={u,4>x}g,..., <f>i= {u,<j>i-i} g , k = l,2, (4.34) 

in order to calculate the desired momenta. Namely, in our case this chain breaks down on the 
third step 03 = 0. It means that H\,i are the second order polynomials in momenta p u and, 
therefore, we can define this unknown momenta in the following way 

Pu = £ (4.35) 

02 

up to the canonical transformations p u —¥ p u + f(u). Similar calculation allows us to determine 
the second momenta p v . 

At k — 1 the results obtained so far can be summarized in the following definition 



Mi 

where 



— (a< - u)(l - du)p. u - (a< - v)(l - dv)p„ . (4.36) 

u — v V / 

(l-<fu)(l-dv) 



(1 — dai)(l — da 2 )(l — dao) 

By adding the expressions for 7$ 



(1-daj )(l-da k ) / u-aQ(v-ai) 
y (1 - rfu)(l - dv) y (aj - ai)(a m - a. ( ) 

obtained from (|3~T8l) to (j4~36"l) one gets the expressions of initial physical variables in terms of 
canonical separated variables. 

By substituting (|4.37[) and (|4.36|) into Hi i2 () 1 . 5[) we can easily prove that variables of 
separation lie on two copies of the hyperelliptic genus 2 curve defined by the following separated 
relation 

4(l-dx)(ai-x)(a 2 -x)(a 3 -x)y 2 -xi? 2 + Hi = 0, x = u,v, y = p„,p„. (4.38) 

It is easy to see, that at d — we obtain the standard elliptic variables on T*S and well-known 
separated relations for the Euler top on the sphere. 
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4.2 Borisov-Mamaev-Fedorov system, k = —1 

Let us take the coordinates of separation qi,2 (|3.27l) and add momenta pi,2 to them, in order 
to obtain a complete set of the canonically conjugated variables (|4.30[) with respect to the first 
Poisson bracket (|2.11l) . 

As above, we can simply calculate these momenta using the recurrence chain <pk (|4.34|) 
associated with the first Poisson brackets {., .} v (|2.1ip at k = — 1 and rewrite the obtained 
results (|4.33|) in the following form 

Mi = ^Moi-c^ (439) 

/ / d(2-b 1 b 2 q 1 ) \ ( d(2-b 1 b 2 q 2 ) \ 
x (Ci - qi) 1 7 Pi - [Ci - q 2 ) 1 7 p 2 



where Cj = ai/bi and 

(l - d(bi + b 2 + b 3 - 2d)qx) (l - d(h + b 2 + b 3 - 2d)q 2 



(1 — dai)(l — dei2)(l — das) aia 2 a 3 



qiq2 ■ 



These expressions are similar to (|4.36p and at d = turn into the standard definitions (|3.8[) of 
the elliptic variables on T 8 S, in contrast with expressions (|4.32[) from [5]. 

By substituting ji (|3.27[) and (|4.39|) into -ffi^ (|1.5p we easily prove that variables of sep- 
aration for the Borisov-Mamaev-Fedorov system lie on two copies of the hyperelliptic genus 2 
curve defined by the following separated relation, 

/_ d(2-b 1 b 2 x) \ (x _ Ci)(x _ C2)(x _ C3)y 2_ ai/2+/3ffi =Q; (44Q) 



63 

where x = qi,2, y = Pi. 2 and 

a = dbib 2 b 3 x 2 - (bib 2 + bib 3 + b 2 b 3 )x + 2 

P = (d( bl b 2 + bl b 3 + W - (bl + &2)(bl + 63)(b2 + &3) ) x + bl + &2 + bs ~ 2d. 

At d = this equation coincides with the separated equation for Hamiltonian systems on the 
sphere separable in elliptic variables, which lie on the elliptic curve instead of hyperelliptic at 
d^O. 



5 Conclusion 

Using the standard machinery of the bi-Hamitlonian geometry, we reproduce some results from 
[5j [9] obtained in framework of the Chaplygin hamiltonization method. Definitions of the 
Poisson bivectors for the Borisov-Mamaev-Fedorov system (I2.11[) . (|3.21D and of the second bivec- 
tor for the Chaplygin case (|3.16l) are completely new. They can be considered as nontrivial 
deformations of the Turiel and Benenti constructions associated only with nonholonomic dy- 
namical systems. 

The explicit form of the separated relations (|4.40[) with canonical variables of separation 
for the Borisov-Mamaev-Fedorov system is also new. Because all the hyperelliptic genus 2 
curves are isomorphic to each other, we could use these separated relations (|4.38l) and (|4.40l) in 
order to get a mapping between Chaplygin and Borisov-Mamaev-Fedorov systems. We suppose 
that such mapping may be extended to the case C 2 ^ 0, that allows us to get solutions of the 
equations (11.4[) in generic case. 

We would like to thank A.V. Borisov for genuine interest and helpful discussions. 
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